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We study a system of two nanomechanical resonators embedded in a dc SQUID. We show that the 
inductively-coupled resonators can be treated as two entangled qubits with states that can be read 
from, or written on by employing the SQUID as a displacement detector or switching additional 
external magnetic fields, respectively. We present a scheme to squeeze the even mode of the state of 
the resonators and consequently reduce the noise in the measurement of the magnetic flux threading 
the SQUID. We finally analyze the effect of dissipation on the squeezing using the quantum master 
equation, and show the qualitatively different behavior for the weak and strong damping regimes. 
Our predictions can be tested using current experimental capabilities. 



I. INTRODUCTION 

In recent years nanoelectromechanical systems-^'^ 
(NEMs), nanoscale mechanical oscillators coupled to 
electronic devices of comparable dimensions, have at- 
tracted substantial research effort. A major motivation 
for this effort is the ability to observe quantum behav- 
ior in a macroscopic system under realizable experimen- 
tal conditional!^. Indeed, NEMs today can be fabri- 
cated with vibrational mode frequencies of 1 MHz-10 
GHz and quality factors in the range of 10'^ — 10^, al- 
lowing the quantum regime to be reached at milli-Kelvin 
temperatures for high frequency oscillators^. Possible 
quantum effects in NEMs under such conditions include 
quantized energy levels, superposition of states, entan- 
glement and squeezing2r— . In addition, NEMs are ap- 
plied to high-sensitive detection of masaiS.^^'^, forcei^ 
and displacementiii, electrometers^^, and also to classical 
memory elements^^'i^. 

Observing or changing the state of NEMs requires 
some type of transducer which couples to them. Op- 
tical coupling^^ can be performed, e.g., by a microwave 
cavity^ but is difficult to integrate in circuits and suffers 
from the diffraction limit and heating of the NEMS. Non- 
optical coupling methods are therefore more common in 
experiments today. With magnetomotive coupling^*^, the 
magnetic force on a thin metallic layer on the NEMS 
is measured. Capacitive coupling can take many forms, 
one of which uses a normal or superconducting single 
electron transistor (SET)^i2i. The NEMS changes the 
island charging energy in the SET and hence the tunnel- 
ing rates, which can be read electronically. Other forms 
of capacitive coupling use Cooper pair boxes^^'^^, flux 
qubits^**, quantum point contacts^^ and quantum dots^^. 

An inductive coupling scheme with a potential for 
displacement precision greater than the standard quan- 
tum limit is obtained by integrating a doubly clamped 
micron-scale beam within a superconducting quantum in- 
terference device (SQUID). In a dc SQUID the motion 
of the resonator changes the area of the SQUID loop 
and hence the magnetic flux and the current through it, 
which is then measured. This system was only recently 



implemented21i2^. A more sophisticated design, where 
the dc SQUID, and hence the resonator, is coupled to a 
charge qubit, was also proposed^. For an rf SQUID it 
was found^S that the change in the magnetic flux due to 
the motion of the beam affects the visibility of Rabi os- 
cillations in the SQUID levels. The detection of discrete 
Fock states in a resonator integrated with an rf SQUID 
was suggested in another work^S. 

Squeezed states, originally introduced in quantum 
optics^, are defined as minimum-uncertainty states with 
less noise in one field quadrature than a coherent stated. 
Several methods to generate squeezing in NEMs were 
suggested. Coupling to a charge qubi1>2.iS as means of gen- 
erating squeezing was proposed, while another work de- 
scribed squeezing by periodic position measurement with 
a weakly-coupled detector^. Squeezing in nanoresonators 
can be applied to decrease the noise in force or dis- 
placement measurements to below the standard quantum 
limit, greatly improving the sensitivity of the device^i^^. 

In this work, we present a scheme to create entangled 
quantum bits (qubits) and squeezing in two nanores- 
onators integrated in a dc SQUID. A previous studj^ 
analyzed a similar system but introduced many approxi- 
mations that are difficult to implement experimentally, 
whereas our present study is closer to an experimen- 
tally realizable system. For instance, we do not over- 
look the generally non-negligible self-inductance of the 
SQUID as done in previous work^i, and we assume mega- 
Hertz frequency nanomechanical oscillators rather than 
giga-Hertz frequency resonators, which are difficult to 
integrate with a SQUID. Lastly, we do not require the 
SQUID to be prepared in a high-|a| coherent state in 
order to have squeezing, as the previous study does^i, 
and require instead a thermal equilibrium state, which 
is easier to accomplish. Finally, we consider different as- 
pects of the system and draw conclusions, e.g., on the 
reading and writing processes, that were not advanced in 
previous literature. 

The paper is organized as follows. In Sec. |n]we present 
the system model and its classical Lagrangian and Hamil- 
tonian formulations. We then proceed to quantize the 
Hamiltonian for the non-dissipative case and derive the 
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FIG. 1. (color online) Schematic of the device we study: Two 
nanomechanical beams oscillating in plane are embedded in 
a dc SQUID with area A when beams are at rest. Uniform 
magnetic field B threads the SQUID, and bias current Ih is as- 
sumed. The two identical Josephson junctions on the SQUID 
have phase drops 7i (i = 1, 2). The currents /s and create 
additional magnetic fields B\ and B-i at the resonators. 



effective Hamiltonian. Next, in Sec. Illll we treat the sys- 
tem as a double qubit and explain how one can read its 
quantum state or write on it. In Sec. llVl we put forward 
a scheme for generating quadrature-squeezed states of 
the nanomechanical beams when dissipation is neglected, 
whereas in Sec. |V]we use a quantum master equation ap- 
proach to test the range of validity of our results in the 
the presence of dissipation. Lastly, we discuss our results 
and present the conclusions in Sec. I VII 



II. SYSTEM MODEL AND HAMILTONIAN 

Our system consists of a dc SQUID, shown schemat- 
ically in Fig. [1] in which each arm includes a Joseph- 
son junction and an integrated doubly clamped beam of 
length li and mass rtii that can oscillate mechanically in 
the plane of the SQUID with an angular frequency uji 
(i = 1,2}. The notation we use is similar to the one 
in Ref. l35l . A uniform magnetic field B is applied per- 
pendicularly to the plane of the loop, and a dc bias cur- 
rent Jb flows through it after splitting to I\ and li in 
the lower and upper arms, respectively. Two current- 
carrying wires with currents 1^ and /4 create additional 
magnetic fields Bi and B2 at the positions of the first 
and second beams, respectively. The beam amplitudes 
are much smaller than the beam-wire distance, allowing 
these fields to be approximated as being spatially uni- 
form. For simplicity, the two Josephson junctions on the 
SQUID arms are taken to be identical and their gauge- 
invariant phase changes are denoted by 7i. The critical 
current and shunting capacitance of each junction are 
taken as Ic and C, respectively, and are used to define the 
characteristic junction energy scales: the Josephson en- 
ergy Ej = hlc/'i^e and the charging energy Ec = e^/2C. 
The plasma frequency ujpi = y/TEjEc /h sets the typical 
time scale for the SQUID dynamics. 

The area of the SQUID loop depends on the center 
of mass positions of the nanomechanical resonators, de- 



noted by Xi and defined as zero when the beam is at 
rest, and positive when it is inside the loop. Since the su- 
perconducting order parameter is single valued, we must 
have 

7i-72-^=27rp, (1) 
$ = - Y.{B + Bi}kX., + L{h - h)/2, (2) 

i 

where $ is the total magnetic flux threading the loop, 
A is the loop area when the beams are at rest, $0 — 
h/2e is the flux quantum, L is the self- inductance of the 
loop, and li is the current in its ith arm. The first term 
in Eq. ^ comes from the external magnetic field, the 
second, responsible for the coupling of the mechanical 
and magnetic degrees of freedom, from the oscillation of 
the beams, and the third from the magnetic flux induced 
by the circulating current in the SQUID. 

The kinetic and potential energies of the system are 
functions of four dimensionless variables defined by 7 = 
(71 + 72)/2, (j) = $/$o and C,: = {B + B,)kX,/<^>o. They 
are 

T = j:{S3^^&) + S^i' + ^^^^ (3) 

U = ii;j[-2cos7cos(7r0) - + Eii-^Y^^^C^ + A) 

i * 
+ f (0-6-6 -0e)^], (4) 

where mechanical dissipation was assumed to be neg- 
ligible, and where we define the screening parameter 
/3l = 2L/c/$o and external fiux (j)^ = BA for the SQUID, 
while the dimensionless magnetic field 

V mi a;i<I>o 

and oscillation frequencies fi^ = uji/ojpi are defined for 
each of the beams. The first term in Eq. (|3]) corresponds 
to the kinetic energy of the beams, while the second and 
third terms to the capacitive energy of the junctions. The 
first term in Eq. ([4]) relates to the Josephson junctions 
energy, while the second term is the washboard poten- 
tial term^s. The third term corresponds to the Lorentz 
force on the beams in the classical equations of motion 
(EOMs), and the fourth term to the beams' elastic poten- 
tial, taken to be harmonic, as nonlinear terms are neg- 
ligible at the amplitudes concerned^^. Lastly, the fifth 
term corresponds to the inductive energy of the SQUID. 

The classical EOMs for the four variables 7, and 
6 are the Euler-Lagrange equations for the system La- 
grangian C = T — U . Before writing the Hamilto- 
nian, we expand the potential in series about a mini- 
mum (0,7,61^2)1 around which the system oscillates. 
Under current experimental conditiona22i2^ such a mini- 
mum exists and the well containing it is deep enough to 
accommodate at least a few states for each of the four 
variables in J7. In the series expansion we neglect terms 
higher than quadratic ones in (j) and 7, as we will later 
find these parameters are "frozen" in the ground state. 
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making the contribution of higher-order terms irrelevant. 
With these approximations in mind and due to the po- 
tential being quadratic, the Hamiltonian H reads 



H = T + U 



EEc 2 



where the coordinates qi are given by 



(6) 



where the interaction coefficient reads 



V 



)1/4(1 



,1/4 



(13) 



III. 



READING AND WRITING QUANTUM 
INFORMATION 



91 = 7-7, 12 = 7r( 
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fe-e,) (7) 



(j = 1,2), and the canonicahy conjugate momenta pi are 
Pi = i^'^ / Ec)qi- In addition. 
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(8) 

where r = cos 7 cos (ttc/)) and s = sin 7 sin (ttc/)) were in- 
troduced. We see that the beam oscillations are coupled 
inductively via the V34 term. This coupling can be used 
to generate squeezed states in the beams as we will show 
below. 

The Hamiltonian is quantized in the standard way by 
converting the coordinates qi and their canonically conju- 
gate momenta pi to operators and postulating the canon- 
ical commutation relation — ihSij. In terms of 
creation and annihilation operators qt and pi are 

1^ - -2iw^m4 - 



(9) 
(10) 



and the quantized Hamiltonian is given by 



(11) 



where uji — LOpi^/Vu- We note the frequencies lo^ and 
W4 are the same as the resonators frequencies, wi and 
UJ2, apart from each having a factor due to the magnetic 
field at the resonator. Thus we see that the Lagrangian 
classical memory variables ^1 and ^2 , in complete analogy 
with memory variables in electronic circuits—, become 
memory quanta in the Hamiltonian. 

Taking the same experimental conditions mentioned 
above and tuning B to make r of order one results in 
fiLJi 3> ksT and hi02 3> ksT. Consequently, the first 
and second harmonic oscillators are "frozen" at their re- 
spective ground states. Moreover, since a;i,a;2 ^ u;3,a;4, 
exciting the nanomechanical oscillators will not budge 
them from their ground state. Removing constant terms, 
we are then left with the effective Hamiltonian 



We now wish to employ this system to create entan- 
gled nanomechanical qubits that can be read from and 
written on. We assume the beams are cooled to a tem- 
perature low enough so as to reduce the equilibrium state 
to the ground state for each of the beams. This is pos- 
sible today, e.g. by laser coolingSS or by coupling to a 
superconducting microwave resonator—, even if the envi- 
ronment of the beams, which includes the SQUID, has a 
higher temperature. 

Perturbation theory in the interaction term in Eq. (fT2|) 
gives first-order energy corrections only when jws — a;4| <C 
V /h. Thus, we will henceforth assume the beams are 
identical. A clearer equivalent form for Eq. (jl2l) with- 
out zeroth order degeneracy is obtained by the canonical 
transformation 



as = -^(03 + ai), 



where a\ corresponds to creation of an even mode quan- 
tum in which both beams oscillate in phase, and a\ to the 
creation of an odd mode quantum, where the beams oscil- 
late in anti-phase. With this transformation the Hamil- 
tonian (IT^ is reduced to 



06 = -^(03 - 04), (14) 



H 



{TiLU-s + V)ala5 + {huj^ — V)alaQ 
+ lv{al + {alr-al^{aln 



(15) 



where we note that the last term is very small relative 
to the first two, since V <^ Tiuj^ with the previously used 
experimental parameters. We see that the even mode is 
decoupled from the odd mode in this Hamiltonian, which 
is thus separable to an even and an odd part. 

Applying time-independent perturbation theory, we 
find the eigenstates of the Hamiltonian to first order in 
V /fiuj'i as 



\nm)' = \nm) + ^Wn{n - l)\n - 2,m) 



- yj{n + l){n + 2)\n + 2, m)} - y/m{m - l)\n, m - 2) 
+ y/{m + l){m + 2)\n,m + 2)}, (16) 

with the first and third terms in the parentheses excluded 
for n < 2 and to < 2, respectively, and with the energy 
levels to the same order being 



H = hwsalas + huj4ala4 + V{a3 + al){a4 + al), (12) 



Enm = (?iw3 + V)n + {hujs - V)m. 



(17) 



4 



It is interesting to note that the interaction results in 
quantum entanglement between the states of the two 
beams as 



1 



V n!m! 
1 



K)"K)'"|oo> 



{al + alrial-alrm. (18) 



V2"+™n!r7i! 



The quantum state of the system is read by measuring 
the magnetic flux threading through the SQUID, which 
is done by a current measurement in the standard way2S,. 
The operator for this observable is 



$ = -(B + Bi)hXi -{B + 82)12X2, 



(19) 



where the constant term BA was omitted. With the 
assumption of Bi = B2, we have 



$ = -V2{B + Bi)hXi{a5 + 4), 



(20) 



where the zero-point fluctuation, the resonator displace- 
m ent uncerta inty at the ground state, is defined as Xi — 

yJh/2miUli+2- 

We thus see, as expected, that the measurement of the 
magnetic flux cannot detect the odd mode, since oscilla- 
tions in this mode do not amount to a change in the area 
of the SQUID loop. We therefore set to read and write 
quantum information only in the n quantum number in 
the state \nm)' . Moreover, we note that for the eigen- 
states of the Hamiltonian we have ($) = 0, implying 
that a better observable would be the standard devia- 
tion (A<I>). This is indeed true, with the values of this 
observable on the eigenstates being 



(A$) = {B + Bi)hXi\ {2 



(21) 



enabling us to measure the value of n. 

Having established the reading process, we now set to 
describe how to write quantum information on this sys- 
tem. It would seem the best way to excite the system is 
via a resonant ac current of frequency uj^ + V /h m. the 
external wires that, according to the Hamiltonian ([T5| . 
will pump the beams to their excited state. However, 
such a current will also pump the beams to even higher 
excited states, since the energy level difference is fixed in 
this system. A better method would be to use constant 
currents in the external wires. The addition to the po- 
tential (U) due to such currents, keeping only first-order 
terms in Bi/B and B2/B, is 

Hi = ^^^{-BiX,[4lS - ^e) + hh] 



-B2X2[hf3L - ^Ic{<l> - cl)e)]} (22) 



where constant terms were omitted and only linear terms 
in Xi were kept, owing to the quadratic terms being 
smaller by several orders of magnitude. 



Since reading can be done only for the even mode, and 
the Hamiltonian (jlSp is separable into an odd and even 
component, we do not consider the odd part in Eq. 
and by choosing also Bi = B2, the even part of Eq. 
is 



Hi, 



4V2ttEjIiXi( 



= /(i?i)(a5+4). 



^Bi(a5 +4) 



(23) 



where constant terms were again omitted. The even part 
of the total Hamiltonian H + Hi is therefore 



{hujs + V)ala5 + ^V{ai + (a^)^) + /(Bi)(a5 -I- 4). 

(24) 

The idea of writing is then the following: We create 
a constant magnetic field Bi, which shifts the harmonic 
potential and then let the system relax to its new ground 
state. We then suddenly revert Bi to zero thereby ob- 
taining an excited state of the original system. The sec- 
ond term in Eq. (I24p only creates admixtures in the new 
ground state and old excited states and will therefore 
be neglected here. The Hamiltonian in the differential 
representation is 

^'&i^^ ^"^^(^^ + + ArV(Si)a:, 

(25) 



1, 



where Ai = h/2mi{uj3 + V/h), and the mass mi does 
not correspond to a real mass and will not appear in the 
final result. Apart from a constant, the Hamiltonian (PS)) 
is 



1 



2toi dx^ ' 2 



-I- -mi(w3 + V/hy{x + AxY, (26) 



Ax 



ArV(Bi) 



(27) 

According to the scheme described above, we wish to 
maximize the probability P„„ = |('0o(^c + Aa;)|-0„(x))p, 
with ipnix) being the nth harmonic oscillator wavefunc- 
tion, to obtain the desired state \n) by tuning Bi and 
with it Ax. Using standard results of the quantum har- 
monic oscillator to write the integral {tpQ^x + Ax)\4'n{x)) 
and then using the generating function of the Hermite 
polynomials to find its value for every n, we find the 
maximum value of P„„ is reached when Ax = 2y^Xi 
with the probability then to measure k phonons after Bi 
is removed given by 



(28) 



which is a Poisson distribution with mean n. The P„„ 
function is plotted in Fig. [2l We see that the probability 
drops sharply for small n and evens out for larger values. 
The limit at n — > cx) is 0. Although the writing process 
does not create a pure number state |n), the standard 
deviation in the number of phonons, by the properties of 
the Poisson distribution, is y^, which is reasonably low. 
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FIG. 2. Graph of the maximum probabiUty of the ground 
state of the shifted harmonic potential having a \n) com- 
ponent in the unshifted potential. We have Poo = 1 and 
Pii = . The values for non-integer n were interpolated by 
n! = r(n+ 1). 



IV. SQUEEZED STATES 

Having shown how to read and write quantum infor- 
mation in this system, we now wish to demonstrate the 
possibility of creating squeezed states. In an effort to 
mimic the Hamiltonian of a degenerate parametric am- 
plifier from quantum optics—, we set the external wires 
magnetic fields to oscillate at double the frequency of the 
even mode, namely, 



Bi = Si.oe2'("^+<^>/'')*-Hc.c, 
B2-B2oe2'("^+<^>/'''* + c.c, 



(29) 



where {V) denotes the time average of V, and the oscil- 
lations of V about this average are small since Bi <C B. 

The non-interacting Hamiltonian is found from Eq. 
(Uni) to be 

Ho = {huj3 + {V))ala5 + (hus - {V})ala6. (30) 

Keeping only terms to first order in Bi/B, the addition 
to the potential (jH) due to the oscillating magnetic field 
in the rotating wave approximation (RWA), valid here 
due to the weak damping, is 



B2)XiX2 + BiXl + B2XI}. (31) 



We write this addition in terms of ladder operators, 
take i?i.o — B2fl to be pure imaginary, note that the 
interaction term from Eq. (1151) does not contribute in the 
RWA and eliminate constant terms to find the interaction 
picture Hamiltonian in the RWA to be 



Hi 



S7r'i(Ai 



E jI-^B\B\^q\ 



{al - (4)2}, (32) 



which is the squeezing Hamiltonian. In writing Eq. (p2)) 
we neglected the time-dependent term {V — {V)){a\a^ — 
agOe), since it can be shown to be negligible relative to 
the squeezing term under the assumed experimental con- 
ditions. It is interesting to note that squeezing of the odd 
mode is not possible using this scheme, even if the sign 
of V is changed in Eq. (P^ . The fundamental reason for 
this is Lenz law, which makes the coefficient of ^1X2 in 
Eq. (|3ip positive and thus excludes terms proportional to 



ag in Eq. (|32|) . The coefficient is positive because moving 
both beams in the positive direction costs energy, since 
both movements decrease the magnetic flux. 

We now consider the effect of squeezing in this system 
and devise means to observe it. We assume dissipation is 
weak, and both beams are initially in the ground state. 
Conforming to standard notation, the squeezing param- 
eter is 

167r(Ai)2£j^25|5i,o|t 



9 



(33) 



which is real, and the squeezing operator reads 

^(.g)=exp[g(a2-(4)2)/2]. (34) 

The time evolution of the 05 operator in the interaction 
picture is given by 



S\g)azS{g) = 05 cosh 5 - ^sinhg. 



(35) 



In the rotating frame the uncertainty in the positions 
and momenta of the beams are 



(AX,) = Ai(l + tanh.g)"i/^ 
(AP,) = i?iA-i(l-tanhg)-i/2, 



(36) 
(37) 



where Eqs. ([7]), ([9]) and ([35l) were used. We see that 
we have limited squeezing to below the standard quan- 
tum limit in the positions and unlimited anti-squeezing 
in the momenta, as the even mode is squeezed, while the 
odd mode is not. In addition, with the product of the 
uncertainties being 



(Ax,)(AP,; 



1 + cosh 2g 



(38) 



the minimum uncertainty is achieved only at i = 0, before 
squeezing takes place. 

It is interesting to write the wavefunction for the beams 
in the x representation in the presence of squeezing. This 
wavefunction can be found by using the relation between 
the two-photon coherent states and the squeezed states^^ 
to write for the squeezed state \g) 

(cosh 505 -f sinhgag)a6|(7) — 0. (39) 



After moving to 03 and 04 operators and then to their 
respective differential representation, Eq. (j39|) is solved 
to find a wavefunction of a jointly Gaussian form, namely 



V'(a;i,X2) = Ci exp(- 



2rx\X2 

(TlCr2 



2(1 -r2) 



(40) 
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where Ci is the normalizing constant, r the correlation 
coefficient and cr^ the standard deviation of Xi. These 
parameters are given by 



_ 1 / l+tanhg ',1/4 
^ Ai\/2¥^l-tanho^ 



tliih^^ ' r = -tanh5, 
CTi = (72 = \/2Ai(l +tanhg)-i/2. (41) 



We see that, as expected, ai = ^/2AXi, where Xi 
IS given m Eq. dSni). The factor of \/2 comes from 
\^p{xl,X2)\'^, rather than ip{xi,X2), being the probabil- 
ity distribution. In addition, we note that the squeezing 
of the even mode increases the influence of the odd mode 
on the correlation, producing a negative r and perfect 
anti-correlation when the squeezing parameter goes to 
infinity. 

Lastly, we consider the effect of the squeezing on the 
measurement of the magnetic flux in the SQUID. Using 
Eq. ((20)) . we find the standard deviation of <I> in the 
rotating frame to be 



(A$) = V2{B + Bi)hXie- 



(42) 



which is fully squeezed, while in the lab frame, with the 
definition uJ-^ = uj^ + (V) /h we have 



(A$) = V2{B + Bi)hXi 
X •\/cosh(25) — sinh(25) cos{2uj3t), 



(43) 



which characteristically oscillates between fully squeezed 
values at i = {-k /u)z)n, corresponding to Eq. (|42|) . and 
fully anti-squeezed values at t = (7r/w3)(n -I- 1/2), where 
n is an integer. We conclude that the squeezing effect is 
measurable, and the squeezing parameter can be found 
from the measurements. 



V. EFFECT OF MECHANICAL DAMPING 

In reality, the damping of the beam oscillations is weak 
but nonzero. With regard to reading and writing quan- 
tum information, this is not a problem, so long as the 
reading or writing is performed within a period much 
shorter than the characteristic decay time. The squeezed 
states, however, are measurably degraded even by very 
weak dissipation, as we show in this section. 

Many models were devised for describing dissipation 
in quantum systems^. We choose here to work with 
the quantum master equation. In the interaction picture 
with the Hamiltonian (I32p , the quantum master equation 
takes the form^^ 



dt 



(44) 

Csp[t)^\g[a^~{o))\p{t)], (45) 
CdisP{t) = -^{ricav + ^){[a\ap{t)] + [p{t)a^,a]} 

- ^ricaviia, atp(t)] + [p{t)a, a^}, (46) 



where for brevity we write a instead of 05 and w instead 
of UJ3. In Eqs. (|44II46I) p(t) is the statistical operator for 



the system, Cs and Cdis are the Liouville operators for 
the squeezing and dissipation, respectively, 7 = lu/Q is 
the damping rate of the even mode, where Q is the beam 
quality factor, and Ucav = [e^^^^^'^ — is the average 
phonon occupation in the even mode. 

The system can be equivalently described by the 
Wigner quasi-probability distribution W{a,a*) instead 
of by the statistical operator p(t), where we omit the 
explicit time dependence in W{a,a*) to make the nota- 
tion concise. The parameter a = xi + 1x2 is a complex 
number that is related to the phase space coordinates via 
a — -^^x + ^p, where x and p are the position and mo- 
mentum coordinates, respectively, and A = y^h/2'miuj as 
before. We convert-2*^! Eq. (|44)) to an equation for the 
Wigner distribution to find 



dW{xi ,X2) 

at 



[5(2; 



1 dx 



X2 



0X2 



+ 



2(90:1^1 dx 



^X2) 



hKa. + i)(g 



+ i^)]Wix,,X2), (47) 



where we note that the original Wigner function W{x,p) 
is related to the one used here by W{a,a*) — 
W{xi,X2) = 2nW{x,p). 

Equation (|47|) is seen to be a special case of the Fokker- 
Planck equation with W^(u) corresponding to the prob- 
ability distribution P(u;i), where u — (xi,X2)- Put in 
this form, the equation can be formally written as 



dWju) 
dt 



-V • [F(u)VF(u)] 



Dr 



V^W{u), (48) 



where F = (-(5+^)a;i, {g-^)x2) and Z^o = 57("cai, + i) 
are the force and diffusion constant, respectively. Due 
to the form of the force in Eq. pS)) . we can use sep- 
aration of variables to break this equation into two 
one-dimensional Fokker-Planck equations with solutions 
Wi{xi) and W2{X2), where W{xi,X2) = Wiixi)W2ix2). 
These solutions are given by (i = 1, 2) 



W,{x,) 



<J^[t) 



1 



exp(- 



(i-^)e-M + ^, 
4 ki hi 



h = 2g + 7, 
fc2 = 7 - 25, 



(49) 

(50) 

(51) 
(52) 



where ki are the decay rates. 

Equations (|49j|52p indicate that a steady-state solu- 
tion always exists for Wi{xi) and is given by Eq. p9)) 
with (Ti(i) = Do/ki. This finite distribution width 
corresponds to a saturation in the squeezing and is dif- 
ferent than without dissipation, when the field quadra- 
ture xi is squeezed without limit. For W2{x2) on the 
other hand, we have a steady-state solution only at the 
strong damping regime, 7 > 2g, and this solution exhibits 
a2{t) = \/Do/k2- When the strong damping condition 
is not satisfied, ^2 is negative, there is no steady state, 
and X2 is anti-squeezed as in the dissipationless case^. 
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but at a slower pace since the leading behavior in [S.Xi is 
e'^9-i/'i)t instead of e^* as in the dissipationless case. 

The knowledge of the Wigner function in Eq. 
enables us to calculate of system properties via the 
relation^ 

{{a-{a^Y}sym)^ j dPaa^{a*YW{a,a*), (53) 

where {-Isym indicates the average of all the permutations 
of the ladder operators, and (Pa — dxidx2- Working 
in the rotating frame, the resulting uncertainties in the 
positions and momenta of the beams read (i = 1, 2) 



(AX,) = V2Ai^-+ai(t)2, (54) 

(AP,) = i=?iAr^y|+^^, (55) 

which reduce to Eqs. p6ll37p when 7 = 0. 

We see that the squeezing in the position coordinates, 
already limited to 1/ of the standard quantum limit 
in the dissipationless case of Eq. is limited here as 

weU to < {AXi)/Xi < 1, where we took ricav — 

due to the previous assumption of huj ^ ksT. The 
momenta uncertainties, in comparison, are anti-squeezed 
only in the weak damping regime, "f < 2g, compared with 
being always anti-squeezed in Eq. ([37]), when there is no 
damping. As with the quadrature field X2, the momenta 
anti-squeezing in the weak damping regime has a slower 
rate relative to the dissipationless case with a leading 
behavior of e'^"'*'/^-'* vs. e^* for the dissipationless case. 
The product of the position and momentum uncertainties 
in Eqs. (|54ll55p gives the minimum uncertainty only at 
t ~ and higher values afterwards. 

As in Sec. IIVI we wish to find here the effect of squeez- 
ing on the measurement of the magnetic flux. Using Eq. 
([20I) , we find the standard deviation in the rotating frame 
to be 

(A$) ^ 2\/2(B + Bi)h\iai{t), (56) 

which is squeezed, though only to a finite extent unlike 
the dissipationless case in Eq. (|42|). where it is fully 
squeezed. In the lab frame we have 

(A$) = 2V2{B + Bi)hXi 
X ^Jcos'^{ujt)al{t) + si-a^ {Lut)al{t), (57) 

which oscillates between squeezed values at t — (7r/a;)n, 
corresponding to Eq. ([55)). and squeezed/anti-squeezed 
values, depending on the damping regime, at t = 
(7r/a;)(n -|- 1/2), where n is an integer. In the dissipa- 
tionless limit of 7 = Eq. ^7} reduces to Eq. (j43]). Eq. 
([57| is plotted in Fig. [3] with normalized units for the no- 
damping, weak-damping and strong-damping cases. We 
conclude again that the squeezing effect is detectible and 
leads to reduced variation in the measured magnetic fiux 
in the SQUID. 




2 4 6 8 

(b) I 
2.5; 




0.5: 



2 4 6 8 10 

FIG. 3. Normalized deviation in the SQUID magnetic flux 
A$(f)/A$(t = 0) as a function of t/Ti, where Ucav = 0, 
Ti = 2-K/ij and g = O.lTi. (a) No damping (7 = 0): Lower 
and upper limits are and 00. (b) Weak damping (7 = 5): 
Lower and upper limits are 1/a/3 and 00. (c) Strong damping 
(7 = 3g): Lower and upper limits are V3/2 and Vs. 



VI. CONCLUSIONS 

In this work we have demonstrated that a system com- 
posed of two nanomechanical resonators embedded in a 
dc SQUID can be used as two units of quantum mem- 
ory — a double qubit — and that only the even mode in 
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these two qubits is readable by the SQUID. We showed 
how the state of the beams can be altered, corresponding 
to writing qubits, and proved the amplitude distribution 
of the number states in the resulting state is Poisson dis- 
tributed. We then proposed a scheme to squeeze the 
even mode of the resonators and thus decrease the noise 
in the SQUID magnetic flux. Taking dissipation into 
account, we found a criterion that separates the weak 
damping regime, where a steady state exists only in one 
field quadrature, from the strong damping one, where 
both field quadratures exhibit steady states. We then 
predicted the form of the fluctuations in the magnetic 
flux in the SQUID, by which squeezing can be observed. 

The approximations and assumptions made during our 
derivations hold well for reasonable experimental values. 
For instance, for two identical 8 MHz resonators of length 
25 fim and quality factor Q = 2-10^, an external magnetic 
field of 10 T, beam temperature of 0.1 mK, SQUID tem- 
perature of 20 mK and other parameter values similar to 
the ones in Refs. [13 and[28l we find the energy level dif- 
ferences in the Hamiltonian (|15p to be much larger than 
both ksT and the level widths. Moreover, for the reading 
process, we note that a typical SQUID with a flux sensi- 
tivity of 10~^^o/^/^^z can detect the smallest differences 
in the flux variation given in Eq. (|2ip for different values 
of n. 

Regarding the squeezing, a major question is whether 
substantial squeezing can be achieved within the decoher- 
ence time for the qubits. The decoherence time for the 
resonators here can be made to be at least 5 /j,si^ii^iii, 
while substituting the parameters above in Eq. (j33p 
gives a characteristic squeezing time of Tsq 2 /is. We 



therefore conclude that substantial squeezing is achiev- 
able within the dephasing time. 

The experimental realization of this system will be 
an important demonstration of macroscopic quantum 
behavior and squeezing in a nanomechanical system. 
In addition it can be used for detecting the position 
of the embedded nanomechanical beams with accuracy 
higher than the standard quantum limit. Stacking such 
SQUIDS in series, with the upper arm of the lower 
SQUID being also the lower arm of the upper one, can 
form a quantum data bus^i^, lead to a multi-mode en- 
tangled stated, and possibly multi-mode squeezing^. 
Another application of this system or a close variant of 
it is that of a quantum gate^ acting on the two qubits 
by means of currents in the external wires. A series of 
such quantum gates can form the basis of a nanomechan- 
ical quantum computer— We leave the development 
of these ideas for future studies. 
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